Entanglement between chiral and deconflnement transitions under strong uniform magnetic 

background field 
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Effects of a uniform (electro-) magnetic background field are investigated in the two-flavor nonlocal 
Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model. Temperature-dependences of the chiral order- 
parameter and the Polyakov-loop are investigated under the magnetic field. In the nonlocal PNJL model with 
four-dimensional momentum dependent distribution function, pseudocritical temperatures of the chiral and de- 
confinement transitions are coincident to each other even if the magnetic field is strong. This means that we can 
check the reliability of the distribution function under the strong magnetic field by comparing with lattice QCD 
data, particularly from the viewpoint of the entanglement between the chiral and deconfinement transitions. 
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It is expected that quantum chromodynamics (QCD) has 
fruitful phase structures at finite temperature (T) and finite 
real chemical potential (/^r) from recent theoretical studies. 
The lattice QCD (LQCD) simulation is a powerful method to 
investigate the QCD thermodynamics at finite T with vanish- 
ing /iR. At finite /iR, the sign problem comes up in LQCD 
simulation and then it is not feasible in the ^J,-R,/T > 1 region; 
for example, see Ref. |1] and references therein. Therefore, 
model approaches are widely used to investigate the phase 
structure at finite /xr. Model approaches, however, have the 
large ambiguities in its foundation; for example, see Ref. [2]. 

Because of the above reasons, we can not obtain reliable 
QCD phase diagram by individually using LQCD simulation 
and model approaches at least in the present day. To over- 
come this difficulty, imaginary chemical potential matching 
approach was proposed in Ref. (jsj]. With this approach, we 
compared model results with LQCD data at finite imaginary 
chemical potential (/ij). After that, we can obtain reliable ex- 
tended model and parameters and then we can clarify which 
interactions are relevant or iiTelevant. 

In this study, we pay attention to a uniform (electro-) mag- 
netic background field as a situation to check the reliability 
of the extended model obtained from the imaginary chemical 
potential matching approach. Moreover, huge magnetic fields 
are expected to be generated in heavy ion collisions 10, Et]; 
then, study of the strong magnetic field attracts much more 
attention. Recently, effects of the magnetic field were ener- 
getically investigated in model calculations. In Ref. |6], it 
is shown that the strength of the magnetic field has an effect 
on the chiral crossover but does not on the deconfinement 
crossover so much in the Polyakov-loop extended Nambu- 
Jona-Lasinio (PNJL) model with smooth regularization. In 
Ref. fh, 

the authors consider the entanglement vertex and the 
eight-quark interaction in the PNJL model with smooth regu- 
larization; then the coincidence between the chiral and decon- 
finement transitions are obtained even if the magnetic field 
is strong. The authors of Ref. [7] remark that the entangle- 
ment vertex is enough to reproduce the coincidence of each 



transition and it is clearly shown that the entanglement be- 
tween the transitions is important in the PNJL model In the 
linear sigma model coupled to quarks and to Polyakov loops, 
that coincidence is also investigated |8]. These studies sug- 
gest that the situation with strong magnetic field provides the 
important information of the model construction, particularly 
for the entanglement between the chiral and deconfinement 
transitions. In a recent LQCD simulation, it was reported that 
coincidence of pseudocritical temperatures of the chiral and 
deconfinement transitions is realized JH. 

In this study, we use the PNJL model H [JMll- There 
are different types of the PNJL model and we use the recent 
PNJL model with the four-dimensional momentum dependent 
distribution function (l^ in this work. The nonlocal nature of 
the PNJL model is natural consequence of QCD nature and 
thus we do not need the cutoff in the momentum integration 
in the nonlocal PNJL model in principle. This is the biggest 
difference between the local and nonlocal PNJL model and it 
is an important advantage if we consider calculations beyond 
the mean field approximation; for example, see Ref. [17] and 
references therein. Moreover, the chiral and deconfinement 
transitions are naturally and strongly entangled in the nonlo- 
cal PNJL model without any other additional interactions and 
vertices; for example, see Ref. lITsl IT9II20I1 . 

The Lagrangian density of the two-flavor nonlocal PNJL 
model is 

C = qiip - mo)q + Ant - U{$[A],§[A];T), (1) 

where q is the two-flavor quark field, mo denotes the current 
quark mass and D" ^ d" + iA" = d" + i5^gAlXa/2 with 
the gauge coupling g and the Gell-Mann matrices Xa- The last 
term U is called the Polyakov-loop effective potential. In this 
study, we treat A'^ as the mean field value as A4 = iA'^ — 
rdiag(0a, 4'bi4'c)- The interaction part is expressed as 



'Cint(a;) = Gsia{x)ja{x), 



ja{x) 



d zC{z)q[x + - Tag x- - 



(2) 
(3) 
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where Fa is -r'a=(o....,3) — (l,i75'f') for the scalar and the 
pseudoscalar interaction, respectively. The function C{z) is 



2 



the distribution function here. In this study, we use the fol- 
lowing distribution function proposed in Ref. ifisll as 



C{p')=C{^lv') = 



(4) 



where — + at finite T with fermion Matsubara 

frequency ujn = (2n + l)7rT. The function as(p^) means the 
running coupling constant. Constants dc and Af are automati- 
cally determined to smoothly connect the function at = 
if we set the actual value of F. It is shown that if we accept 
the bare /i-dependence of the distribution function by using 
following replacement of po. 



Po 



iAA 



(5) 



the A4-dependence of C is a natural consequence from ex- 
istence of the Roberge-Weise periodicity [21] as shown in 
Ref. ||20|[ . The A4-dependence is also shown from QCD based 
analysis 1^^. In this study, we neglect the Z factor that comes 
from the quark wave function renormalization to make our 
discussion simple. Details of the Z factor in the present non- 
local PNJL model are shown in Ref. [19., 2Q1. 

In this study, we use the following logarithmic Polyakov- 
loop effective potential lfl5|[ as 



U = F'' 



+ 64(T) ln{l - 6<Z'<P + 4(<?3 + <Z>3) - 2,{<P<Pf} 

&2 = ao + ai(^— j +a2(^— j +<^3[y) ' 

, /Ton 3 

04 a4(^Y. 



(6) 
(7) 
(8) 



where the <P and <P with the Polyakov-gauge are <P = ^e**^" + 
e#i, + e**=)/3, ^ = ^* here = (03 + 08/V3)/T, 0^ = 

(-03 + 08/\/3)/T and 0, = -(</)„ + 0b). 

From here, we consider the uniform magnetic background 
field A^^*. The magnetic field along the positive z axis is 
expressed as AJ^^* = Bxi5^2 where we choice the Landau 
gauge. In this case, the four-dimensional momentum becomes 
P= (PO,0,Q v/2fc|Q/eS|,p^) and then p2 ^p2_^2fc|g/eS| 
below where k labels each Landau level and Q denotes the 
sign function as sgn((5/) with electric-charge Qf for each 
quark flavor f — u,d. Making the mean field approximation 
and introducing Landau levels, the thermodynamical potential 
can be expressed as 
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where I3^l/F,E)^ ^f- + (M})2, cj„,, = (2n + l)7rr + 
{(^{F - ifi) and M = diag(M„, Md) with = = 



mo — 2GsCi{p^)(T here i represents the color indices. It is 
shown that the /i is set to zero in this calculation. The co- 
efficient I3k = 2 — Sko express the degeneracy of the Lan- 
dau level. The thermodynamical potential (|9]l is difficult to 
converge with numerical summation over n. Therefore, we 
rewrite it as 
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where SIq means the thermodynamical potential with nonin- 
teracting quarks but it is affected by the uniform background 
field Aa and A"^^* . The final form is 



no 



dpz 
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ln/ + ln/ 



(11) 

/ = l + 3{^ + 3e-'^'''^P'>)e-'^'''^f'> +e-^'^^'^P\ (12) 
/ = 1 + 3(<P + (Z.e-'''=^(p))e-'^'=^(P) + e-3'3«^(P\ (13) 



with = y^fp+m^zt /i. The vacuum part of J?o is already 
subtracted from thermodynamical potential (fill because it 
does not affect the position of the global minimum. 

The present nonlocal PNJL model has three parameters in 
the NJL-part, Gs, mo and F. The Polyakov-loop effective 
potential has six parameters, oq, oi, 02, 03, 04 and Fq. Pa- 
rameters in the Polyakov-loop effective potential are fixed by 
using LQCD data in the pure-gauge limit. The actual values 
are uq = 3.51, ai — —2.56, a2 = 15.2, 03 = —0.62 and 
04 — —1.68. The energy scale Tq is set to 270 MeV. Parame- 
ters in the NJL-part are fixed by using empirical values of the 
pion mass and its decay constant. Actual values are mp = 3.3 
MeV, Gs = 20.6 GeV-^ and F ^ 0.827 GeV. 

First, we show results of the PNJL model with smooth reg- 
ularization to compare with the nonlocal PNJL model with 
a four-dimensional distribution function. Here, we do not 
present actual formalism of it; see Ref. [7] for the details of 
the formalism and the parameter set. Below, we call the PNJL 
model with smooth regularization and the four-dimensional 
distribution function as 3D-PNJL model and 4D-PNJL model, 
respectively. Figure [T| shows the F dependence of the chiral 
condensate and the Polyakov loop in the 3D-PNJL model for 
eB/jnl = 0, 10 and 15. In the case of the 3D-PNJL model, 
the Polyakov loop is less sensitive than the chiral condensate 
against the strength of the magnetic field. This means that 
the difference between pseudocritical temperatures becomes 
large when the magnetic field becomes strong. Rough estima- 
tion of pseudocritical temperatures are summarized in Table |I] 
In the estimation, we use the temperature-derivative of order 
parameters, and the numerical error is about ±0.001 GeV. 

Next, we show results of the 4D-PNJL model. Figure |2] 
show the <T field and the Polyakov loop for eB/m^ = 0, 10 
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Fig. 1 . The T dependence of the chiral condensate and the Polyakov- 
loop in the 3D-PNJL model for eB/m^ = 0, 10 and 15, respec- 
tively. The chiral condensate is normalized by that at T = with 
eB = 0. 
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Td (4d-PNJL) 


0.208 


0.215 
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TABLE I. Rough estimation of pseudo-critical temperatures for the 
chiral (T^) and the deconfinement transition (Td) determined from 
the peak position of the temperature derivative of order parameters. 
Unit of and Td is GeV and these have about ±0.001 GeV numer- 
ical error 



ulation l'9'l. The 4D-PNJL model can naturally reproduce the 
LQCD prediction without some other additional interactions. 
At /i = with eB = 0, transitions are considered as the 
crossover or the very weak first-order by recent LQCD data; 
then it is difficult to clearly check the reliability of the entan- 
glement effect in the model. However, the transitions become 
sharp first order under a strong magnetic field and then it is 
possible to clearly check the reliability of the model through 
the dependence of the model against the strength of the uni- 
form magnetic background field. 

Finally, we discuss our results with the data in Ref. JtII- 
Figure |3] shows the percentage of increase of the pseudocriti- 
cal and critical temperatures as a function of eB/rn^. In the 
figure, Tc means the pseudocritical or critical temperatures 
for each eB/m^ where shows the pseudocritical temper- 
ature with eB/nri^ = 0. The solid line represents the results 
of the EPNJL model picked up from Fig. 6 in Ref. |01 and 
this model coiTesponds to the 3D-PNJL model with the entan- 
glement vertex. We only plot the pseudocritical temperature 
of the chiral transition because pseudo-critical temperatures 
of the chiral and deconfinement transitions are very close to 
each other in the EPNJL model. Opened and closed sym- 
bols are the pseudocritical and critical temperatures of the 4D- 
PNJL model, respectively. The 4D-PNJL and 3D-PNJL mod- 
els show good agreement and the crossover is changed into 
the first-order above eB /ml. ^ 11 in the present 4D-PNJL 
model. 

1.2| ■ ■ ■ ■ , ■ ■ ■ ■ 1 



and 15. The a field can be considered the chiral order pa- 
rameter which is the indicator of the chiral transition with fi- 
nite quark current mass because it appears in the quark mass 
function. In the case of the 4D-PNJL model, the pseudocrit- 
ical temperatures of the chiral and deconfinement transitions 
are very close to each other at several eB / - the same as 
te LQCD suggestion and actual values are also shown in Ta- 
ble|T] Moreover, we can see that the transitions are close to the 
first-order behavior when the magnetic field becomes strong. 
This first-order tendency is already suggested in LQCD sim- 



Our result 
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Fig. 2. The T dependence of the chiral order parameter (cr) and 
the Polyakov-loop in the 4d-PNJL model for eBjm^ — 0, 10 and 
15, respectively. The chiral order-parameter is normalized by that at 
T = with eB = 0. 



Fig. 3. The percentage of increase of the pseudocritical and critical 
temperatures as a function of eB/m\. Here, Tc means the pseudo- 
critical or critical temperatures for each eB /m^ and T^ means the 
pseudocritical temperature with eB/rnl- — 0. The solid line repre- 
sents the EPNJL model results shown in Fig. 6 of Ref. |7]. Opened 
and closed symbols are the pseudocritical and critical temperatures 
of 4D-PNJL model, respectively. 

In the 4D-PNJL model, the form (|5]) is the natural (minimal) 
form by respecting the replacement manner of pq, the preser- 
vation of properties of QCD at finite imaginary chemical po- 
tential [21] and manifestation of the recent QCD based anal- 
ysis ilSll when we accept the bare ^ dependence in the dis- 
tribution function. Importance of this form is also expressed 
in the functional renormalization group analysis ||22| . There- 
fore, there is the possibility that we can indirectly investigate 
the ^ dependence of the model even at /i = through Eq. (|5]l. 
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Of course, the form (|5) is a sufficient but not necessary con- 
dition. Therefore, we must carefully investigate the /i depen- 
dence step by step. 

We can use the system with the strong uniform magnetic 
background field as the situation to quantitatively check 
the model reliability in the near future when more accurate 
LQCD data are obtained. 
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